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Skyrmion gases and lattices in helimagnets are studied, and the size of a Skyrmion in various 
phases is estimated. For isolated Skyrmions, the long distance tail is related to the magnetization 
correlation functions and exhibits power-law decay if the phase spontaneously breaks a continuous 
symmetry, but decays exponentially otherwise. The size of a Skyrmion is found to depend on a 
number of length scales. These length scales are related to the strength of Dzyaloshinskii-Moriya 
(DM) interaction, the thermal correlation lengths, and the strength of the external magnetic field. 
An Abrikosov lattice of Skyrmions is found to exist near the helimagnetic phase boundary, and the 
core-to-core distance is estimated. 
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I. INTRODUCTION 

Skyrmions are non-trivial topological objects in var- 
ious field theories. They were first used to model 
baryons in nuclear physics.'^ More recently, they have 
been observed in quantum Hall ferromagnets,^^ p-w ave 
superconductors,'* and Bose-Einstein condensates.'^ The 
Skyrmion lattice is also a candidate for an experimentally 
observed periodic phase in helimagnets such as MnSP^^ 
(called the A phase) and Fei_2,Coa;Si (called the SkX 
phase). *° The schematic phase diagram of a helimagnet 
is shown in Fig. [l] In these helimagnets, the usual or- 
dered pha se tha t exhibits long-range order takes the form 
of a helix,'iiJi2! ^j^^t is thermodynamically stabilized by 
Dzyaloshinskii-Moriya (DM) interaction.*'^ The helix 
of the magnet can be characterized by the pitch vector 
q, with magnitude being the helical wavenumber and the 
direction being the direction of the helix. Its magnitude 
is proportional to the strength of the DM interaction. 
The phase that is believed to be a lattice of Skyrmions 
in these helimagnets have lattice size of the order of mag- 
nitude q~^. It is likely that dilute Skyrmion gases in var- 
ious phases of helimagnets can also be stabilized through 
the DM interaction. 

An isolated Skyrmion is an exact solution of the saddle- 
point equation to the non-linear a model.'^-^ *'^ Thermody- 
namically it is a metastable state. Different schemes have 
been proposed to stabilize the Skyrmions, with a dipole- 
dipole interactiorP^ and a DM interactiorP^JtSil being two 
examples. These schemes introduce new physical length 
scales. 

A lattice of Skyrmions has also been describe d as 
the superposition of three perpendicular helimagnets .1^*251 
This Fourier description was argued to be a lattice of 
Skyrmion by calculating the Skyrmion density. However, 
this approximation is not a solution of the saddle-point 
equation to the model they are using.E^ Using the tech- 
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FIG. 1. Schematic phase diagram of MnSi in the H-T plane 
showing the helical, conical, and A phases, as well as the 
field-polarized ferromagnetic (FM) and paramagnetic (PM) 
states.liS 



nique of Abrikosov vortices in type-II superconductors ,^21 
a lattice of Skyrmions can be described. The Goldstone 
modes for such lattice is the same as that of the colu mnar 
phase of the liquid crystals."^ Some previous studie^^SHU 
showed that by increasing the magnetic field, the system 
changes from aligned conical phase, then a Skyrmion lat- 
tice and finally to a ferromagnet. This is true for certain 
temperatures but not for all. The Skyrmion lattice is 
likely formed by a first-order phase transition. 

The sizes of Skyrmions can be characterized in four 
ways. It can be characterized by the behavior near 
the core.l^ Alternatively, the size can also be character- 
ized by the decay length of the long-distance tail of the 
Skyrmion. In this paper, we introduce yet another mea- 
sure of the size of a of Skyrmion, which is given by the 
length at which the behaviors of the core and the tail 
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match. We believe that it is the best measure of the size 
of a Skyrmion. Finally, in a Skyrmion lattice, there is a 
fourth measure of the size of a Skyrmion: the distance 
between cores in the lattice. 

In this paper, we study both an isolated Skyrmion and 
a Skyrmion lattice in the presence of DM interaction and 
external magnetic field, and we estimate their sizes. In 
Sec. In) we review the Landau-Ginzburg- Wilson (LGW) 
model with the DM interaction used to study the heli- 
magnets. Crucially, we also introduce the various length 
scales in the problem. In Sec. |III[ we review the ba- 
sic properties of Skyrmions, including the winding num- 
ber. In Sec. |IV) we study the core of isolated Skyrmions 
in various phases of the model. We find that for para- 
magnets and ferromagnets, the core behavior is readily 
found. The core behavior defines a measure of the core 
size, called R, which we find to be the result of the com- 
petition of different physical length scales in different re- 
gions of the phase diagram. We show that for the aligned 
conical phase, the core behavior is undetermined because 
of the scale invariance due to the spontaneous symmetry 
breaking. In Sec. [Vj we study the Skyrmion tails in var- 
ious parts of the phase diagram. For paramagnets and 
ferromagnets, we find that the tail is exponentially de- 
caying. The decay length of the tail depends on various 
length scales in different regions of the phase diagram. 
The Skyrmions in aligned conical phase, however, are 
algebraically decaying for large distances due to the un- 
derlying Goldstone modes in the system. In Sec. |VI[ we 
employ the Abrikosov flux lattice to study the Skyrmion 
lattice. We find that the core size depends on q and the 
thermal correlation length of paramagnets. In Sec. |VII[ 
we match the core and tail of Skyrmions, with the match- 
ing radius introduced as a new measure of the Skyrmion 
size. We argue that this measure of Skyrmion size is the 
most physical. In general, we find that the core size is 
of the order near the helimagnetic phase boundary, 
consistent with the results of other previous studies, ^'^ 
but it depends on other length scales in other parts of 
the phase diagram. We do not discuss pinning, polariza- 
tion, and the alignment effects mentioned in a previous 
paper.li^ 

II. MODEL AND LENGTH SCALES 

A. LGW functional 

Throughout this paper, we use the LGW model with 
DM Interaction^^ 



5[M] = 



^(VM)2 + |m-VxM 



-(M2)2 - H-M 



(1) 



The terms with coefficients r, a and u are the usual 
Landau-Ginzburg- Wi lson (LGW) model.'^Sl The c-term 
is the DM interactio d^^ l '^'^l , and it exists in systems with 
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FIG. 2. Phase diagram predicted by the action in Eq. ([T|. 

spin-orbit couplin^SHl and no inversion symmetry in the 
unit cell of the solid.'^ A similar chiral structure can be 
found in cholestoric liquid crystal.'^i' H is the external 
magnetic field. Its direction defines the z-axis through- 
out the paper. The saddle-point equation corresponding 



to Eq. (jlj) iS' 

rM - aV^M + cV x M + uAf^M - H = 0, (2) 

We do not consider crystal-field effects in this paper. 
Generally these terms align the helix in a certain direc- 
tion for small magnetic field. For large magnetic fields, 
the helix is in the direction of the field. 



B. Phases and Length Scales 

There are three topologically trivial phases that are as- 
sociated with the action in Eq. Q. These three phases 
are the paramagnet, the ferromagnet and the aligned con- 
ical phase, as shown in Fig. ([2]). The aligned conical 
phase is stable only if _ff < Hc2 where Hc2 is defined 
below, see Eq. ( 18 ). 



1. Paramagnet 

The paramagnet is stable or met a st ab le only for r > 0. 
Its analytic form is given by Eq. (B2). In a magnetic 
field, the magnetization is 



XpH, 



(3) 



where Xp — is the magnetic susceptibility. The corre- 
sponding correlation length is given by ^p, which diverges 
when r approaches 0. In LGW model, it is^^ 



(4) 



We define a length (basically the thermal and magnetic 
field dependent transverse correlation length) 



(5) 
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For r > it is 



(6) 



But for r <C to « {u ^H)^. Hence the length be- 
comes 



r— T, 



(7) 



which is related to the mean-field critical exponent 6 = 



2. Ferromagnet 

The ferromagnet is stable or metastable only for r < 0. 
Its analytic form is given by Eq. ( B4 ) . Without external 



magnetic field, the magnetization is given by 



(8) 



A magnetic field along the direction of the magnet gives 



M 



AO) 



XfH, 



(9) 



where Xf = (2k|)^^ is the longitudinal susceptibility. 
The translational susceptibility goes like H~^. 

The longitudinal correlation length of the ferromagnet 
is given by 







(10) 



which diverges for r — )• 0. The transverse correlation 
length is given by Ih in Eq. ([s]). It is infinite for zero 
magnetic field, but in a field it is 



H 



(11) 



where Xh = {a.q^)~^ is the magnetic susceptibility of 
the homogeneous magnetization, nisp is the helimagnetic 
amplitude given by 



TOc 



where 



niH = 



aq'^ — r 



(15) 



(16) 



The aligned conical phase can be the thermodynami- 
cally stable phase as in Fig. ([2|, if the system is in the 
region 



2 uH^ 
aq -r> 

This defines the critical field Hc2 



Hc2 = aq' 



2 aq-r 



(17) 



(18) 



At zero field, the correlation length approaching the 
helimagnetic-paramagnetic phase transition from the 
paramagnetic phase (r > aq^) is given by 



1 



r - aq^ V^p 



q 



(19) 



If the phase transition is approaching from the helimag- 
netic phase (r < aq^), it is given by 



(20) 



At the transition point between the helimagnet and para- 
magnet at = 0, 

qip = 1. 



III. BASIC PROPERTIES OF SKYRMIONS 



3. Aligned Corneal Phase 



A. Winding Number 



The aligned conical phase, or simply conical phase, 
contains a component of helicity and another of homoge- 
neous magnetization, given by^^ 

Mcp — TOsp(cosg'^x + singzy) + to//z, (12) 

where q is the pitch vector of the helix given by 

It is proportional to the strength of DM interaction, m// 
is the homogeneous magnetization induced by the mag- 
netic field, and is given by 



V/ = XhH, 



(14) 



Skyrmions are two-dimensional objects, which are 
topologically non-trivial because the winding number of 
a Skyrmionic configuration is non-zero. Assume that 



M(x) = m(x)n(x), 



(21) 



where to(x) and n(x) denotes the magnitude and direc- 
tion of M. The winding number is defined aJ^ 



W ~ I dx I dy^n ■ 



dn dn 

dx dy 



(22) 



Upon continuous deformation of the configurations, the 
winding number W remains unchanged. Note that all 
phases in Sec. |II B| have = 0, which means they are 
all topologically trivial. 
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B. Description of Skyrmions 



(23) 



To illustrate a topologically non-trivial solution, we 
write the configuration in the form o^^ (in cylindrical 
coordinates) 

n(x) = sin0(x) cosa(x)/5 + sin0(x) sina(x)<y9 + cos0(x)z, so that the winding number can be written as 



1 f'OO p27T 

— dp d(p -smOix) 
47r Jo Jo 



96l(x) aa(x) 
dip dp 



9a(x)\ 56* (x 



dip J dp 



(24) 



For all the Skyrmions we review and present in this pa- 
per, 9{p = 0) = TT and 9{p = oo) = 0, so that W = —1. 
In addition, we set a = ^, meaning it is an azimuthal 
Skyrmion. We exclude the consideration of a radial 
Skyrmion because the presence of DM interaction forces 
the Skyrmions be azimuthal. The winding number de- 
scribes how the electron changes its spin when it passes 
through the core.^ Such a Skyrmion is depicted in Fig. 
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3. The picture of an isolated azimuthal Skyrmion with 
-1. 



Although Skyrmions are topological objects, linear re- 
sponse ensures that the asymptotic behaviors of the tails 
of Skyrmions are no different from the decay of other 
kinds of perturbations. It is the characteristics of the 
core that constitutes the topology. In understanding 
Skyrmions in various ordered phases of the helimagnets, 
we study the core by differential equations with bound- 
ary conditions that give a non-zero winding number, and 
the tail through perturbation theory. 



0) = TT. Then we write, for p — > 0, 

m(x) = moo + 6ra{p), 
0{x) =Tr + Se{p), 



(25) 
(26) 



where all behavior depends only on the radial distance 
p, and where TOqo = Mp (for paramagnet, r > 0, as 



in Eq. ( |B2[ )) or Mp (for ferromagnet, r < 0, as in Eq. 
(|B4l)). Putting Eq. ([25| and Eq. ^ in the saddle-point 



equations (Dl I and (D2), ignoring all higher order terms, 
we can find that for p ^ 0, 69 — Q^P and 6ra — M^p'^, 
with 



= - 



moo 
1 + 4q% ' 



(27) 
(28) 



where Ih is the transverse correlation length for param- 
agnets or ferromagnets in Eq. (|6]), Eq. ([7) or Eq. (Ill 



Note that this calculation breaks down if there is no DM 
interaction, or the Skyrmion is in an ordered phase. 

The size of the core can be estimated with this solution 
as R, where 9 = 7r(l— ■^). Let us consider different 
cases in ferromagnets and paramagnets, and the relation 
of Skyrmion in paramagnets in some situations to the 
Skyrmion lattice. 



1. Ferromagnet 

For a ferromagnet with non-zero magnetic field, the 
core size R is given by the following cases (with the phase 
stable or metastable): 



1. Ih < 



R 



2q- 



IV. SKYRMION CORES 

A. Ferromagnet and Paramagnet 

We first explore the cases of paramagnets and ferro- 
magnets. And as we have stated in Section III 9{p — 



2. < Ih: R^2'Kql'\j. 

The first case refers to the region where the magnetic 
field is much larger than the critical field Hc2- The core 
size is proportional to q~^ . The second case refers to the 
region closer to H = Hc2, where the magnetic length Ih 
plays a role. 
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2. Paramagnet 

For a paramagnet, the core size R is given by (in which 
the phase can be stable or metastable) : 



(D5| and (D6|, the dominant variation near the core is 



1. £.p < 



R. 



2q- 



2- ~ 9 ^ and along H w Hc2'- R ■ 



Stt 
2q- 



The magnetic field does not play a role in the core be- 
havior for paramagnets. It is because close to r « 0, the 
paramagnet appears for H > Hc2 and Ih is less signifi- 
cant than the contribution of For r > Ih ~ S,p- 



B. Aligned Conical Phase 

The core behavior for aligned conical phase is different 
from paramagnets and ferromagnets in Sec. |IV A| because 
of the Goldstone mode in this phase. To understand 
it, we first study ferromagnet without external magnetic 
field and DM interaction. This is for illustrative purpose, 
because it breaks a continuous symmetry just like the 
conical phase does. 



psin^. However, the coefficients is undetermined for the 
same reason as Sec. |IVB 1| The perturbation for conical 
phase can be written as Eq. (41). The differential equa- 



tions for the fluctuations for conical phase can be written 
as a Laplace equation as in Eq. ( [42| , although both gra- 
dient terms and curl term (due to DM interaction) in the 
action in Eq. ([T]) are important for the aligned conical 
phase. 



V. SKYRMION TAILS 

The tail behavior of the Skyrmions is closely related 
to the correlations in the bulk of the phase. This can 
be shown by linear response theory (see Appendix |e]) . In 
the following, we show by each case that the presence of 
spontaneous symmetry breaking gives the Skyrmion an 
algebraic tail, but an exponentially decaying tail other- 
wise. The decay length of the tail is in general not the 
same as the core size. 



A. Paramagnet and Ferromagnet 



1. Ferromagnet without external magnetic field and DM 
interaction 

With H — and q — 0, the ferromagnet (with r < 0) 
is an ordered phase that breaks the continuous rotational 
symmetry of the action in Eq. ([T]). The differential equa- 
tion for 69 in Eq. ( 26 1 as in Eq. (|D2| becomes 



dp'^ p dp 



(29) 



The differential equation is due to the gradient term in 
the action. Then 



= Ap, 



(30) 



with an arbitrary coefficients A. In fact, the approxima- 
tion in Eq. (C7) from the exact Skyrmion solution to the 



non-linear a model near the core gives the same. There- 
fore, we cannot determine the core size. To determine 
the core size just from the core behavior, external per- 
turbations that provide extra length scales are needed. 
For exam ple, the presence of DM interaction fixes I in 
Eq. ([C7| to be (4g)-i. 



We study the system using perturbation techniques. 
As we have stated in Section III 6{p 
we write, for p 



oo) = 0. Then 



■> oo, 

m(x) = 
6i(x) =04 



) + Sm{p), 

seip), 



(31) 
(32) 



where all behavior depends only on the radial distance p, 
and moo = Mp (for paramagnet, r > 0, as in Eg . ( |B2[ )) 
or Mp (for ferromagnet, r < 0, as in Eq. (B4)) as in 
Sec. |IV| We expect that the tail is exponential. Hence, 
we assume Sm = Me~^P and S6 = Qe~^P. One of the 
solutions for K, called K^, is used to define the length 
of the tail 



1 



(33) 



which is another measure of the size of Skyrmions. In 
some cases, has an imaginary part, which indicates 
the tails are oscillating in addition to the exponential 
decay, but we will omit oscillations below despite its ex- 
istence in some cases. 



Ferromagnet 



2. Aligned Conical Phase 

The aligned conical phase breaks the continuous trans- 
lational symmetry of the action in Eq. ([l]) Because of 
the helical nature of this phase, we expe ct m (x), ^(x) 
and q;(x) depend on ^ (defined in Eq. (D3)) in addi- 



tion to the radial distance p. After analyzing Ec^s. (D4|, of ferromagnet have the spectrum w(k) = k 



For ferromagnet, the lengths of the Skyrmion tails are: 

1. Ih < < if. h ~ Ih- 

2. Ih < < q^^ and ^/ < < q^^: h ~ Ih- 

For both cases, the tails has a length It ^ as in Eq. 
(|ll|). We expect this because the transverse fluctuations 



2. Paramagnet 



For paramagnet, the lengths of the tails are: 



1. Ih < < ^p: h 



H- 



3- ^ and small i/: « 



4. Along w i/c2: 



where for the fourth case, 



K± 



\ 




(34) 



The first case corresponds to the the paramagnet with 
large magnetic field and the boundary with ferromagnet, 
where Ih is given by Eq. The second case refers 

to the paramagnet far away from the transition points, 
making It ~ £,p- The third case refers to the paramagnet 
very close to the helimagnetic transition point, making 
It - a- 

The fourth case refers to the paramagnet along H w 
Hc2- The Skyrmion lattice is formed along part of this 
critical line. 



From an analysis of the differential equations (Dl| and 



(Dl I, we get 



0/ 

P ' 



Sm 



,r(0)^ 



(37) 
(38) 



From Eq. (36), we know that the coefficients 8/ is ar 



bitary, as in core solution the coefficients in Eq. ( 30 1 in 
Sec. IIVB H is undetermined as well. It is the result of the 
equation for the fluctuations given by a Laplace equation 
for a symmetry-breaking ferromagnet. 5m{p) and 59{p) 

can be expressed in the scaling form / {^jj, £^ , for some 

length scale i. Moreover, the approximation in Eq. (C7) 



far from the Skyrmion core in the non-linear a model has 
the same behavior. Similarly, as in the Skyrmion core, 
the coefficients can be fixed by additional interactions 
that carry other length scales as discussed in Sec. |IVB 1| 



2. Aligned Conical Phase 

While the ferromagnet has the readily diagonalized 
Goldstone modes with spectrum w(k) = fc^, the aligned 
conical phase has the Goldstone modes given bj^ 



B. Aligned Conical Phase 

By Goldstone theorem and linear response, any pertur- 
bation in the conical phase shows long distance algebraic 
decay. Therefore, a Skyrmion in this phase shows a long 
tail. For illustrative purpose, we study the ferromagnet 
without magnetic field and DM interaction which breaks 
the rotational symmetry. 



1. Ferromagnet without external magnetic field and DM 
interaction 

The spectrum of the Goldstone modes in ferromagnet 
in = and q — Q \s a;(k) = fc^,'22! and the modes are 
readily diagonalized as Snix and 5my. As a result, they 
behaves |r|~^, as illustrated in Appendix [Ej We expect 
the Skyrmion tail to behave in the same way. 

We still employ the perturbation schemes in Eq. (31 1 
and Eq. ( 32 ) for this ferromagnet with H = and q = 0. 
The differential equations kept to the relevant order is 
given as 



q2 



AO) 



H-2 



H 

Hc2 



H 

c2 



(39) 

for small magnetic field {"rn^/ / ^ ''^h)^ similar to that of 
the cholesteric liquid crystal.!^ For zero magnetic field 
makes, the second term vanishes.'^ For the parametriza- 
tion of fluctuations about the conical phase ( 12 ) can be 
written aJ^^ 

M(x) = TOsp [cos{qz + ',5o(x))x + sin{qz + <po(x))y 
-f((^+(x) cos qz + </?-(x) sin(7z)z] 
-|-m//[7ri(x)x + 7r2(x)y -I- z], (40) 

to the first order of all parameters. The Goldstone mode 



that corresponds to (39) is given by the "diagonalized" 
form as 

(5mG(x) 







9(^+(x) 


q{q'^ 




dy 



H 
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which satisfies the partial differential equation 



H 

Hc2 



5mG(x) « 0, (42) 



for H ^ 0.-^^ Eq. (42) has a solution 



SmcM = — sin((/3 + B), 
P 



(43) 



with undetermined coefficients A and B. (Note that the 
Skyrmion core in conical phase has arbitrary size R and 
goes like pshi(ip + B) for the same argument.) There- 
fore, similar to ferromagnets, the Skyrmion tail in the 
conical phase has a power law form. A detailed analy- 
sis of the saddle-point equations in Appendix |D 2| shows 
that Skyrmion tail in the conical phase goes like p^^ sin^, 
where ^ is defined in Eq. (D3 1, with arbitrary coefficients. 
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Hexagonal Skyrmion lattice described by (45 I, for 



on the plane, and the color denotes Uz — z*al-' Zj where deep 
blue denotes spin reversed from the magnetic field. 



VI. SKYRMION LATTICE 



In this section, we study a lattice of Skyrmion in he- 
limagnets. It is convenient to use CP^ representation 
because Skyrmions are like vortices in superconductors 
and the technique of Abrikosov lattice of vortices can be 
used. The details of this representation can be found in 
Appendix |A] The saddle-point equations in this repre- 
sentation are Eq. ( |A6| and Eq. ( A7 1 . 

In Eq. (A7), the gauge A depends on z as in Eq. 
(A4), leading to the non-linearity. To fix the gauge, set 
A — —hxy. The dimension of h is that of the reciprocal 
of area, and it will be shown later that it is related to the 
area of a Skyrmion core. Because of the periodic nature 
of the lattice, the term —iqmriadaZi — —iqm{na)daZi is 
ignored. This can be justified by the final solution. The 
term —^qmzidaria is also zero because the Skyrmions are 
azimuthal. Moreover, instead of keeping strictly z'^z = 1, 
we relax the condition to (z^^z) = 1 where the average is 
over one lattice. Following AbrikosoVjI^Sl one part of the 
solution is given bj^^H 



1 



h 

1 

+ Mol 



1 

2h (x 



where the prefactor is for normalization, and 
iq 1 



(44) 



do 



2h 1 

2 



H 

Aahm 



1 4. 31 
4 ^ 2h 



H 
iahm 



' V Aahrn I 



Let Ix and ly be the distances between cores along the 
X and y axes respectively, where Ixly = . Then z can 
be seen as the superposition of the above solution with 
different values of k where fc,- = ^ , then 



h 
34 7r 



j = -oo 



For a triangular lattice, Cj 



1 

do\/2h [x + jlx) 
(45) 



Choose Co 



for even and odd fs. And ly 



2 



to be 
Such 

2 



configuration is plotted as shown in Fig. |4] for h ^ q 
which denotes a Skyrmion lattice. On the other hand, a 
graph with spin projected on the basal plane (x-y plane) 
and a density plot of = z*a^J Zj is plotted in Fig. [ij 

To know the magnetization and the core size, we have 
to put Eq. ( 45 1 back to Eq. (|A5|) and determine them 



by variational method. Since the lattice is periodic, it is 
valid and convenient to consider the solution of a single 
site in Eq. (|44|. Define 



h 



V 




the free energy per unit volume of one Skyrmion in the 
lattice is given by 



V 



h 



2 2 

-aq m - 



1 



2h 



2V 1 _ 
Hm 



V 
2h 



1 



2h 

(46) 

Then we evaluate magnetization m and the reciprocal of 
core area h by minimizing the free energy. There exists 
no analytic closed form solution for m and h, but we do it 
by qualitative analysis. We expect that m is of the same 
order of magnitude of m// in the aligned conical phase 
or the paramagnet, and ^ < 1. Expanding Eq. (46) for 
small we get 

V 2 



— Hm 



2 2 

aq m 



Aa?q^ 
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Minimizing it with respect to h and m, we get 

h i/2 



rm + uni^ — H fa 0. 

The second equation indicates that the magnetization is 
approximately equal to the paramagnet or ferromagnet. 
For small magnetic field and in the paramagnetic phase 
(r > 0), using Eq. (B2), we get 



1 



2t4 



(47) 



Near 



Therefore, the core-to-core distance goes like q£_p. 
the helimagnetic phase boundary ^pq « 1, and in this 
region the core-to-core distance goes like q~^. 

The experimentally observed A phas^ that has been 
identified as hexagonal lattice of Skyrmions is obse rved 
along the helimagnet /paramagnet phase boundarypESl 
Our result shows that the core size ~ q^^, which is ~ q~^ 
for g^p « 1, agreeing with prev ious theoreticaP^^MIllIIIl 
and experimental studies-l^ESl However, we also predict 
that the size increases away from the phase boundary, 
provided the Skyrmion lattice is still the thermodynamic 
ground state when the correlation length gets larger. 

The fluctuations due to the Skyrmion lattice has the 
same form of that of columnar phase in liquid crystal,'^^! 
and have the form a;(k) = k^^ -I- cifc^ + C2kj^ 



VII. CORE SIZE AS THE MATCHING 
DISTANCE BETWEEN CORE BEHAVIOR AND 
SKYRMION TAIL 



The definition of R in Sec. |IV A| deals only with the 
core behavior. Here we introduce a new distance L, which 
is defined as tbe distance where the core behavior and tail 
of the Skyrmion meet. By finding L, we consider both the 
core and the tail of the Skyrmion. In some cases, R and 
L are not too different in terms of order of magnitudes, 
but their difference becomes greater when the magnetic 
field becomes large. 

In the following cases, we match the core behavior and 
the tail at a point L, and then we solve for L. 



A. Paramagnet and Ferromagnet 

For paramagnet and ferromagnet, by matching the 
core behavior in Sec. |IV| and the tail in Sec. |V] we 
solve for the matching point L tabulated in Table^ We 
verify that for all cases in Table |l] have winding number 
W = —1, by putting the solutions of 0(x) back to Eq. 



(24). 



From Table |T| we can see that i?'s are mostly of the or- 
der of magnitude of q~^, and Zt's are mostly the thermal 
correlation lengths. However, L shows much more com- 
plicated dependence on the various length scales. It is 



generally not a good measure of a Skyrmion size because 
the correlation length is related to the thermodynamic 
phase of the bulk, the size of an additional object. Both 
R and L is of the order of magnitude of q~^ near the 
helimagnetic transition point at iJ « 0, indicating that 
the Skyrmion size is of q^^ in this region. Far from this 
point, L and R differs in orders of magnitude. In general, 
L is better to characterize the size of Skyrmions because 
it takes into account both the core and the tail. Whether 
L or i? is a better measure depends on the situations, as 
listed below: 



1, 



2. 



If L and R are of the same order of magnitude (~ 
K^^ or oscillating) as in Fig. [5](a), they are equaUy 
good. Examples: paramagnets in q^^ Ih £,p 
and < q^^. 

If i? L as in Fig. [s] (b) , L is a better measure be- 
cause L depicts where the tail starts and the slope 
of the core behavior was underestimated. Exam- 
ples: ferromagnets in £_f <^ q~^ <^ In , £.f ^ Ih ^ 
Ih < < and Ih < q^^ < and 



paramagnets in Ih 'ti q ^ 
Ih- 



< Cp and ^p < q-^ < 



3. If L ^ i? as in Fig. [5] (c), i? is a better measure. 
The matching method is not working so well be- 
cause at p ~ L, 6{p) becomes negative. However, 
There are no such examples in all cases considered 
in Table m 

Therefore, for our purpose, L is a better characterization 
of the Skyrmion size in general. 




FIG. 5. Plots for 9{p) for Skyrmions in paramagnets and 
ferromagnets, and the meanings of R and L as the sizes of the 
Skyrmion core, (a) R fa L. Both are equally good measures 
for the core size, (b) R ^ L. L is a better measure for the 
core size, (c) R <gi L. 7? is a better measure for the core size. 

For pure ferromagnets {q — Q and H = 0) and aligned 
conical phase, matching does not fix the size of the 
Skyrmion due to the same reason stated in Sec. |VB 1| 
For pure ferromagnets, matching the solutions in Eq. 
(C7) for small and large p does not give L. 
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TABLE I. Isolated Skyrmions in ferromagnets and paramagnets in different regions of the phase diagram in Fig. [2] (PM = 
paramagnet, FM = ferromagnet) 



PM/FM 


Region in Phase Diagram 


It 


R 


L 


FM 




Ih 


2%qFu 


^■i2nqH%if 


FM 


C/ < Ih « q'^ 


Ih 


2q 


/ T ^ 

^JA-KqlHt,f 


FM 


Ih < « q'^ 


Ih 


2q 


\J2-KqlHlH 


FM 


I'M y ^-ts. <~,j 


Jh 


TT 

2q 




PM 


< q'^ < Ih 


Cp 


TT 

2^ 


1 


PM 


q'^ < < Cp 


-27-1 

q Ih 




2TvqpH 


PM 


Ih < q'^ < Cp 


Ih 


TT 

2q 


yj2TtqlHlH 


PM 


Cp < q-' 


a 


57r 
2q 




PM 


Along // « _H'c2 


V2 


Stt 
2q 


TT 



B. Aligned Conical Phase 



where (in terms of bare parameters in the action) 



Because of the technical complexities of aligned conical 
phase, we do not match the core and tail solutions as we 
did for other magnets. However, we assume a form of 
Syrmionic solution here and estimate the range of the 
core size. We describe an isolated Skyrmion in terms of 

M = [1 - rii{p,iy\map{cosqzyL + sinqzy) 

+mi/nsk - V3{p, 0'7^//z, (48) 



where risk is the same Skyrmion in (Eqs. [Clp3| . The 
first and third term are with variational parameters rji 
and 773 which does not alter the winding number. The 
second term with hsk keeps the winding number to be 
— 1. Because Eq. (48) already captures the long-range 
behavior as hgk ■ (— sing2;x + cosgzy) ^ ^ sin^ far from 
the core (p — > 00), there is no term in the direction of 
(— sin qzx + cos qzy) in Eq. ( 48 ) . 



We then solve for rji and 773 by putting Eq. ( 48 1 to the 
saddle-point equation ([2| for regions far from the core 
(p 00) and near the core (p w 0). Far from the core 
(p — >■ 00), we find that 



El 



[r{l-lq) + l§{2-q^P)] m 



II 



r Y- + u{Smlp + 4m^^) 

2w(™// - 2TO2^)m//TO2^ 
(7m2^-4m2p) r - + u{3mlp + 4m'j ^) 
uql{28rrt'il - 23m^/^mlp + 4:m//m'^p) 



(53) 



{1m?i I — 4m2p) r — ^ -|- w(3m^p -|- 4to^^) 



S3 = 



\qlrr? 



sp 

Amlp 



Aa 

;2 



^ + r - Aurf^ii + um\p 



2a 



4TO2p' 



(54) 
(55) 



{Zaq^ - 2r)r 



V/ 



21 



sin^ 



m/ / (21 



irisp \ p 



sin^ , 



(49) 



(50) 



where I is an undetermined parameter that we estimate 
below. At the core {p « 0), with linearization of the 
parameters in the saddle-point equation ([2]), we find that 



2p 



m~-Ei[^smC] , (51) 



(52) 



Then we match Eq. ( 49 ) and Eq. ( |51[ ) , and 773 by match- 
ing Eq. (50) and Eq. (52) at some distance p — L, and 
we can solve for I. Matching is only possible if B3 and E3 
are positive, and Ei is negative. For to be positive. 



the denominator Imj ^ 



iml has to be positive. (This 



makes the Skyrmion gas to appear only if > 0.798i?c2-) 
L depends on ^ slightly. The lack of an analytic solution 
forces us to explore the core size numerically. 

In Table |ll) the sizes of Skyrmions in the aligned con- 
ical phase with different values of the correlation length 
£_'i^ and H = 0.8iJc2 are listed. L is the size of a Skyrmion 
found by matching method. It depends on the phase an- 
gle ^ but it does not vary significantly. In general, as 
the system goes away from the phase boundary (as 
decreases), the size of the Skyrmion decreases. 
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TABLE II. Isolated Skyrmions in aligned conical phase for 
different values of all for H — 0.8Hc2- 





I 


L 


£ — 02 — (D 
S y 


3.16g-^ {r = O.Yaq'^) 




Z.ooq 


U 




Z.oiq 


TT 

2^ 


1 .oyy 


z.ooy 


TT 


1 c;q^— 1 
i .Ooy 


Z.Z4y 


Stt 
2 


1.83g'^ {r = 0.7oq^) 






U 


1 . lOy 


1 .ooy 


TT 

2 


1.13g"^ 


2.13g"^ 


TV 


1.12g"^ 


1.38g~^ 


Stt 
2 


1.41g~^ [r = 0.5ag^) 


0.88g~^ 


1.35g~^ 





0.88g~^ 


1.29g"^ 


TT 

2" 


0.88g-^ 


1.35<j"^ 


TT 


0.88g-^ 


1.23g"^ 


Sir 
2 


1.29q-^ [r = 0.4ag^) 




1.13q-^ 







1.15g-i 


TT 

2" 


0.79g-^ 


1.13g-^ 


TT 


0.80g~^ 


1.20<3~^ 


Stt 
2 



It is also of order q ^ near the helimagnetic phase bound- 
ary, and decreases away from it. 
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Appendix A: Model in CP^ Representation 

To facilitate the study of topologically non-trivial 
Skyrmions, it is useful to write the model in terms of 
the CP^ representation, which is commonly used in the 
study of quantum Hall ferromagnets.!^ We first write 



M(x) = m(x)n(x). 



(Al) 



where m(x) and n(x) denotes the magnitude and direc- 
tion of M. The direction n can be written in terms of 
two-component spin through the Hopf mappingPSl 



(A2) 



VIII. SUMMARY AND CONCLUSION 

In summary, we have studied various aspects of 
Skyrmions in hehmagnets. We first studied the Skyrmion 
core, and a related measure of core size R for para- 
magnets and ferromagnets. For the aligned conical 
phase, this R becomes arbitrary. Extra perturbations are 
needed to determine it. We then discussed the Skyrmion 
tails. They decay exponentially for paramagnets and fer- 
romagnets. We determined their decay lengths It, which 
is another measure of the size of Skyrmion. We can- 
not define It for the aligned conical phase because the 
Skyrmions have algebraic tails, as expected from Gold- 
stone theorem. We also studied the lattice of Skyrmions. 
Through variational methods, we found that the core- 
to-core distance in the lattice is of the order of magni- 
tude q^p. This can be compared with the experimen- 
tally observed A phase where the core-to-core distance 
is found to be ~ q~^. Insofar as q£,p « 1 near the heli- 
magnetic phase boundary, experimental and theoretical 
results agree. Lastly, we introduced the matching radius 
L, as another measure of the core size. In our studies, 
the core sizes are of the order q^^ near the helimagnetic 
phase boundary, consistent with various previous stud- 
ies. However, in other parts of the phase diagram, it 
depends on other length scales as well such as the ther- 
mal correlation lengths or and the magnetic length 
Ih , as shown in Table |l| Among all measures of sizes of 
Skyrmions, we think that L is the best because it de- 
pends on the whole Skyrmion. We also estimated the 
size of Skyrmions in aligned conical phase by matching. 



where cTq is the Pauli matrix. From now onwards, Greek 
indices denote spatial component and Latin indices de- 
note merely matrix component in spins. The constraint 
= 1 gives 



Zi Zi = 1. 



(A3) 



Because n has a definite magnitude, there are only two 
degrees of freedom. Therefore, z should have two de grees 
of freedom only, accomplished by the constraint Eq. ( A3 ) 
and fixing the gauge^^ 



(A4) 



such that the transformation ^^(x) — > e^'^(''^Zi(x) and 
z*(x) — 7> e'^^^'^z* (x) does not lead to a change in the 
physical system. 

In this representation, we write the action in Eq. ^ 

as 



S[m, A, z] 



= / d\ 



2 4 2 

+2am^ [{daZ*){daZi) - AaAa 

1 



damdaiTi — Hamz*al^ Zj 



-cm 



z*a'^ZjAa, + '-{do^zDa'iz, - '-z*a'Jdc.zi](^) 



The saddle-point equation associated with this represen- 
tation has to be done wi th L agrangian multiplier because 
of the constraint in Eq. (A3 1. Consider 5 -I- A / (Px{z*Zi — 



1) with A being the Lagrangian multiplier, the saddle- 
point equations are 
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(r — 3aq^ — 4:aAaAa)m — adadam — AaqmUaAa + um^ — Haria+Aam 



iq 



= 0, 
(A6) 



TO (^ijda - iSijAa + ^CTa ^ (^SjkOa ~ iSjkAa + ^'^a''^ ~ iqm (^UadaZi + ^Z^daUa^ + ^^aCr'a = Az^. (A7) 



Appendix B: Analytic Expressions for the 
Paramagnet and the Ferromagnet in Mean-field 
LGW Model 

For the paramagnet and the ferromagnet, the magneti- 
zation can be solved by the saddle-point equation derived 
from LGW functional, 

rm -I- mn^ - H ^0. (Bl) 

Its analytic solutions are given in the following. 



1. Paramagnet 

The paramagnet phase is only valid for r > 0. Its full 
analytic expression is given by 



Mp = XpHf 



(B2) 



where 



1 

M /x\ 2 ^ 1 X /l /x\ 



(B3) 

A Taylor's expansion for small H confirms that Mp w 
XpH as in Eq. (|3|. However, when r approaches 0, Alp « 

(^) ^ , giving the critical exponent (5 = 3 for the mean- 
field theoryP 



2. Ferromagnet 

The ferromagnet is only for r < 0. In general the 
magnetization is given by 



M 



(0) 



Mp = 



(cos- 



^/3^ 



where Mp^ is given by Eq. ([8|), f{x) is defined in (B3) 
and 



27\r\ Y kl 

m 



sin7 : 



27 nr 



2\r\ 



rH. 



(B5) 



A Taylor's expansion for small H confirms that to 



M 



(0) 



Appendix C: Isolated Skyrmion in Non-Linear a 
Model 



In the ferromagnetic phase, the fiuctuations can be 
studied with non-linear a model.l^ In this model, there 
exists a metastable mean-field solution that corresponds 
to a Skyrmion given bjff^ 



2Zpsin(^ 

" P + p2 ' 

2lpcos if 



p2 + /2 ■ 



(CI) 

(C2) 
(C3) 



where I is an arbitrary length that characterizes the size 
of th e S kyrmion core. In terms of the representation in 
Eq. (pi, 4 



6l(x) = 2 tan- 
In CP^ representation. 



p^ + P 



P' 

I 
p 



(C4) 
(C5) 

(C6) 



which is an anti-vortex.'22l This solution has a winding 
number W = —1,^^ confirming that it is a Skyrmion. 
This Skyrmion in non-linear a model states that 



y -, for p ^ oo. 

We cannot determine I as discussed in Sec. IIV B II 



Appendix D: Saddle-Point Equations for m(x), 6(x) 
and a(x) 

1. In Field-Polarized Magnets 



Putting M = mn with n given by Eq. (23) with 
a(x) — I into Eq. ([2|, the saddle-point equation be- 
comes 



XfH. 



/d^m ldm\ / 89 
\ dp^ P dp J \dp 



1 . 2 

77 m sm 

P^ 



12 



' d0 1 

-cm I — — \ — sin 6 cos d \ + unC' — H cos = 0, 
.op P 



(Dl) 



139 1 ^ ^ 2 dm 39 

am + - t; 5- sm cos ^ — — 

\3p'^ pop p'^ m op Op 

+c — V -m sin^ 9] - iJ sin 6* = 0. 

\dp P J 

(D2) 

These two equations are useful for describing a Skyrniion 
in paramagnets and ferromagnets. 



2. In helimagnetic phases 



Putting M = mn with n given by Ec^. ( 23 1 and with 



the Eq. ^ gives the following equations: 



(D3) 



rm, sm a cos a ~ a 



2 3 



ld_ 
p3p 



p—ims sin 9 cos a) 
3p 



TT^im^ sin 9 cos a) H — ^ sin sin a) ?r?7is sm9 cos a 



- ^ [ms cos 9) - q^{ms sin 6^ sin a) 
p3^ 3S, 



umi sin 9 cos a = 0, 



(D4) 



rm, sm f sm a — a 



1 9 



1 



p3p 
{ms sin 9 cos a) 



p—ims sin 6* sin a) 

9/3 



TO, smfsma 



— — (m„ cost^) + q — lms smt^cosa) 
3p^ ' 3C 



+ um sin 9 sin a = 0, 



(D5) 



13_ 
p3p 



+c 



Pi^inis cos 9) 
3p 

1 9 



rTO^, cos 9 — 
13 

- — ipms sin 9 sin a) -\ — — (m^ sin 9 cos a) 



3s;- 



:{ms COS 9) 



+ um „ cos t 



H = 0. 



(D6) 



Appendix E: Linear Response Theory 

Linear response theory relates external perturbations 
of a physical system to the correlation functions in the 
unperturbed one!^ This ensures that the perturbations 
have long-ranged behaviors if there exist correlation func- 
tions for massless modes for fluctuations. In contrast, 
the perturbations are of short-ranged if there do not ex- 
ist any massless modes. This explains the behaviors of 
Skyrmion tails in different phases in Sec. |V] in relation 
to the Goldstone theorem. 

Assume there is a field M(x) with a known mean- field 
solution Mq, with the corresponding free energy per tem- 
perature being Sq. Let the kernel matrix for fluctuations 
be K(x, x') and an external perturbation H(x). Then 
the corresponding partition function is 



Z = e 



-So 



DM- 



exp 



■ J d''^ J d'^x'^SMi{x)K,j{x,x')6Mj{y:') 



(El) 



where summation is on repeated indices and /3 = 
{ksT)^^. By carrying out the functional integral, we 
get the fluctuation determinant .1^^ The equation for 6M, 
by variational principle, is 



and hence, 



)nce, 

(5M,(x) = 13 J • K-:\x,x')Hj{x.'). 



(E2) 



(E3) 



Note that the inverse of the kernel matrix is actually the 
correlation matrix. We can see from here that the pertur- 
bation (5M(x) has the same behavior as the correlation 
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functions K(x,x'). If ii'-Hx,x') = ivf-^x - x'), the This gives SM^{x.), SMy{^) 



Fourier representation of Eq. ( E3 ) is 



and SM;,{-x.) 



(E4) 



We illustrate this with the example of a pure ferromag- 
net with q = and H = 0. Assume that the magnet is 
aligned along z direction. Then the kernel matrices (the 
inverse of the correlation functions) are given bjEH 

Kij{k) = Sij{aP - 2rS-iz). 

This indicates that along x and y directions, the fluctu- 
ations are massless as expected by Goldstone theorem. 
Along z direction, it is massive. Then, we get 

ak^SMyik) - PHy{k) =0, 
(ak^ - 2r)5M^{k) ~ /3iJ^(k) = 0. 



|x|~^e . Here, we have demonstrated that the cor- 
relation function of a massless mode leads to long-range 
behaviors of the perturbations. Similar behavior exists in 
the perturbations to aligned conical phase in Sec. |VB 2| 
Here, wc illustrate an example of a pure paramagnet 
with q = and H = here. In this case, Mo = 0. The 
kernel matrices ar^^H 



K,,ik) 



Sij (ak^ 



Then 



[ak^ - r)8M,{k) - /3i?,(k) = 0, 



for i — x,y and z. This gives 5]V[(x) ^ |x|^^e '^p . Here 
we have demonstrated that if the correlation functions of 
all modes are massive, then the perturbations are short- 
ranged. 
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